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ABSTRACT 


We report our prepress ia computing theoretical solar chromospheric and 
photospheric models for use in analyzing OSO 8 spectra. We have updated the Vemazza, 
Avrett, and Looser (1976) solar model and produced self-consistent non-LTE number 
densities for H I, He I, He H, C I, Mg I, A1 1, Si I, and H". These number densities 
have been used in the calculation of a theoretical solar spectrum from 90 to 250 nm, 
including approximately '7000 lines in non-LTE. ‘More fhan 60, 000 lines of other 
elements were treated with approximate source functions. 
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TIEE3 CALCUIATION OF 

THEORETICAL CHROMOSPHERIC MODELS AND 
PREDICTED OSO I SPECTRA 

NASA Grant NSG-7054 
Progress Report No» 1 


1. INTRODUCTION 

The OSO 8 e^iperiment has observed radiation from different regions on the solar 
disk and, through specstrmn scans, has probed a range of depths in the atmosphere 
from the corona to the upper photosphere. Those observations, together with earlier 
Skylab and rocket observations, should lead to a new and greatly improved understand- 
ing of the structure and behavior of the solar atmosphere and of the heating of the 
corona, as well as to significant advances in our theoretical understanding of high- 
temperature plasmas. 

The variation with depth of temperature, density, and other atmospheric param- 
eters at a given point on the disk can be determined from the emergent spectrum. 
However, the computations necessary to derive these parameters recpiire a major 
effort. We have developed a two-stage procedure for computing theoretical solar 
chromospheric and photospherlc spectra. In the first, we construct a solar model 
based on the non-LTE computer program PANDORA (see Avrett and Looser, 1969; 
Avrett, 1971; Vemazza, Avrett, and Looser, 1973). 

In the second stage of the spectrum calculation, we make use of a spectrum- 
synthesis program written by Kuruez (1974) and based upon his model atmosphere 
program ATLAS (1070), which includes data for 300, 000 lines computed by Kuruez and 
Pesrtremann (1975), The statistical-equilibrium number densities from the model 
calculation are used to determine the non-LTE line source functions of the abundant 
elements, while the less abundant ones are treated with approsimate source functions. 
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la tho noxfc seoticm we describe our current solar model,* and in section 3 we out- 
line various developments in the spectrum synthesis program, especially the treatment 
of lines merging into continua and non-LTE modifications to the basic programs. Ih 
sections 4 and 5 our spectrum calculations are described, finally, we discuss further 
work to remedy shortcomings in the most recent calculations. 
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MODEL 


We construct a solar model by oboostog a temperature-height distribution such 
that the computed emergent intensity is in best agreement with observed line and con- 
tinuum spectra in the ultraviolet, visible, and infrared spectral regions. PANDORA 
is used to solve the radiative-transfer and statistical-equilibrium equations ’lydrogen, 
helium, carbon, silicon, and other important constituents. The model calcui £ion is 
iterative in that only one atom can be computed at a time, while others are held fixed. 
Thus a self-consistent solution is found by cycling through the atoms several times, 
updating the number densities with each calculation. 

Vemazza, Avrett, and Looser (1976) have completed a detailed stucfy of tlie solar 
continuous spectrum redward of 125 nm. They have constructed a temperature-density 
model, which is a representation of the quiet solar atmosphere in the height range 
between the photosphere and the chromosphere-corona transition region. We have 
revised and improved that model to obtain ilie temperature-height relation shown in 
Figure 1. The model is self-consistent in H I, He I, He H, C I, Mg I, A1 1, Si I, and 
H". 


We ejq)ect to fiirfcher revise this model by studying the differences between our 
calculated spectrum and the rocket and OSO observations all through the ultraviolet 
and the Slqrlab spectra below 140 nm. 
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PiguTG 1. The temperature-heiglifc lelEtioa for our revised solur model* 


3. SPECTRUM SYNTHESIS PROGRAM 


We have conttaually improved our spectrum synthesis programs during the course 
of this research. The basic spectrum synthesis program works as follows. For a 
given model atmosphere tabiilated at 40 depths, we compute number densities and parti- 
tion functions for all atoms and ions. Then we s:ead the magnetic tape that lists the gf 
values calculated by Kuruoz and Peytremann (X976), and compute a line opacity spec- 
trum at each depth. The expressions for radiative, Stark, and van der Waals broaden- 
ing are 


^R “ ^classical 223E13A^ , X wavelength in nm 5 

rg=l.0E-Sn^/N^ i 

= 4 . 5E-9 Mg + 0.42 Njjg + 0. 86 NHj! (T/10000)®'^^® j 


where 







i 


is the charge plrs X, E^^p is the upper level energy, and R is the Rydberg energy. 
Next the continuum opacities, source functions, and intensities are evaluated at several 
points in the v^velength interval. The continuum quantities are interpolated to each 
point in the spectrumj the line opacity is added; the source fimction is approximated 
hy 

[(/c + + i (K + CT + i ) ; 
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and the emorgont Bpoctrom is computed. Here k and a are tto continuous aijsorption 
and scattering ooeffioients, i is the line absorption coefficient, is ttio continuum 
source function including scattering, aM is tlio lino source function, wbicli is 
assumed to bo the Plonolc function. 

The following improvements have been made. We have computed mean square 
i?adil for the 4p orbitals of the iron group atoms and ions, and produced an interpola- 
tion formula 

<r^>= (45-A + Z)/(Z+l) 

whore A is the atomic number and Z is the charge. This e:q>resslon is used for 
van der Waals broadening of all iron group atoms, assuming that the strong linos are 
all resonance lines with a 4p upper state. 

There is now provision for additions, deletions, and corrections to the table of 
gf values and for use of individual "real" radiative, Stark, and van der Waals damping 
constants for each line. 

To resolve the contribution of individual lines to blonds, we now compute, in 
addition to the spectrum, the central intensity for each lino in isolation. 

We automatically produce spectrum plots with identifications. 

The two most important advances as far as improving the appearance of the cal- 
culated spectrum have been treatment of lines merging into continua and the develop- 
ment of a procedure for computing thousands of lines in non-LTE. 

One of the obvious faults in the preliminary calculations made one year ago was 
that strong discontinuities appeared at the 01, Si I, A1 1, and Mg I edges, whereas 
the real spectrum makes smooth transitions from one continuum to another because 
in the real speetjum there are numerous lines that gradually merge into each contiu- 
uum. In rewriting the opacity routines to include departures from LTE as discussed 
helow, photc ionization continua were put in for level-to-level transitions rather than 
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for term-to-torm transitions as ms tho case previously. TIius, for example, for C I 

q O 

and Si I, ionization from the neutral ground term P to the ionized ground term P 
actually forma Oeoutaa, 

^p 2 ^^2 "*’*^3/2* those continua has a series of lines converging 

up to it. The lines merge into a quasi-continuum at an upper energy or upper quantum 
number determined by Sfark broadening, because lines are no longer visible when tlio 
Stark width exceeds the lino spacing. Jn practice, at each depth in the atmosphere, 
we compute the energy at which individual lines are no longer visible and e^end each 
continuum redward at a constant level to this point. This technique smooths the 
transition from one continuum to another even if no linos are included explicitly. 

Por A1 1 and Mg I we were able to compute Rydberg series to high n (80), and to 

find or estimate the radiative. Stark, and van der Waals broadening for each line. 

For Si I and C I there is so much configuration interaotiem that extrapolated energy 

-1 

levels can be off by up to 2 cm imless great effort is made to include the Interaction 
effects, so we were not able to compute the series. For Si however, energy levels 
have been observed up to n » 66 although many of them were unclassified. We were 
able to pick out the interacting levels by estimating the mutual repulsion, so v/e could 
assign each level to an LS series. Then we estimated LS f values by extrapolating 
each series. These f values will be correct for many of the lines but wrong for the 
strongly interacting lines. However, the total f values spread over all transitions 
should be correct as configuration interation only redistributes the line strengths. 

To demonstrate the improvement in the new calculationa we show in Figure 2 
sample old and new spectra calculated around the Si I edge at 198 nm. Tlie old cal- 
culation predicted an abrupt discontinuity, and the new one, a gradual slope, in agree- 
ment with the observations. 

To treat the non-LTE effects we have rewritten ATLAS (Kuruoz, 1970) routines 
for partition functions and for continuous opacity from H I, He I, He H, Cl, Mg I, 

A1 1, and Si I to produce non-LTE populations, opacities, and source functions. 

We have added a new line opacity routine that produces nrai-LTE opacity and a source 
function for each C I, Mg I, A1 1, and Si I line. The changes were accomplished by 
inserting departure coefficient b*s as follows. 
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Figure 2* (a) A previous calculation showing the silicon discontinuify at 197* 8 nni* Also shown are two Sirvard 
rocket sc^ by Kohl, ParMnson, and Seeves that do not show an abrupt disccmtinuity. The vertical 
scale is luiear but in arbitraiy units* Our most recent calculation ^cswing tbe smooth traiisItioB. 
across this region, in agreement with the observations, which results from treatment of lines mersinu 
into the continuum. 




partition functions: 




-Bj/kT 

0 


wlioro and aro tho statistical weight and onorgj^ for each level: 
Bolt 2 anan foctors; 



where N and Nj are tho total and level number densitiosj 
Stimulated emission: 

J.g-hV/lCT_J.Vg-hvAT 

”lo 


where up f io refer to the upper and lower levels. For ionization, up refers to the 
ground state of the ion. 

Source functions; 


S(lo^up) = 






b. 


lo hv/kT 

r"® 

up 



In practice, b's are calculated only for the lowest 6 to 10 terms, so for higher levels 
b’s are assumed to have ihe same value as the b for the lowest ls%'el of the ion. 
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4. COMPUTING TIME AT NOAB 


This work on spectrum synthesis has been exceedingly difficult, requiring inno- 
vations and much more sophisticated programing than originally envisioned. Further- 
more, calculation of the spectrum in non-LTE requires a great deal of computer time, 
even for debugging. Because of limited funds provided by the grant for computing, 
we have to use free computing where available. Fortunately we have been able to use 
the CDC 7600 at the National Center for Atmosplierio Besearoh throu^ collaboration 
with OranB. White of the High Altitude Observatory. A considerable inconvenience of 
working at NCAB is that to check the calctilations we must bring magnetic tapes back 
to Cambridge to plot the spectra because NCAB does not have plotting facilities. 
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5. spectrum CALCULATION 


Wo used the model descriljQd in section 2 and the program described in section 3 
to compute the theoretical solar spectrum from 90 to 2S0 nm including non-LTE effects 
for H I, He I, C I, Mg I, Ai I, Si I, and H". We computed approximately 7000 lines 
in non-LTE. More than 60000 lines of other elements were treated with approximate 
source functions. 

In Figure 3 we show an enlargement of Figure 2b with line identifications. Of the 
many lines included in the calculation, only those with residual intensities less than 
0. 99 are indicated. There are numerous non-LTE Si I lines in this region. The lines 
that show up prraninently wiih triangular shapes in both the observed and calculated 
spectra are Si I resonance lines. 

In Figure 4 we show a small region of chromospheric spectrum at 149, 3 nm where 
the lines appear in emission. 


U 



Figure 3* The calculated spectrum compared td tlie Harvard rocket spectra of Kohl, 
Parkinson, and Reeves for the region 197. 6 to 199. 1 nm. The observed 
and computed spectra are not plotted on the same vertical scale; the two 
scales are linear hut arbitrasy. This figure shows approximately the same 
TOvelength region as Figure 2b but enlarged to shew detail and with line 
identifications added. The identifications at the top of the figure include; 
the last three digits of the wavelength, the element identification and charge 
(atomic number before decimal point, charge after: thus 14. 00 = Si I and 
28. 01 = Ni E), then, the lower energy level in cm“^ and finally the fraction 
of the Ml scale in per mil for the center of each line if it "were in isolation. 
Thus strong lines have small central intensities and weaii lines have values 
near the continuum, which is indicated. The numerous Si 1 lines are com- 
puted in non-LTE. 
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Figure 3. 
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Hgure 3 (Continued). 
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Figure 3 (Continued). 
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Figure 3 (Continued). 
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6, 2?t3ETHER WORK 


We expect to malto a numlber of iinprovemonts over the next two months, then 
return to NCAR for another iteration on the spectrum. Several problems that we 
expect to correct are -the following. 

C H, Mg H, A1 n, and SI H lines are computed too weak because they wore not 
based on non-LTE results, Work is underway to calculate departure coefficients for 
these ions, 

MaJTJf highly ionized lines do not appear in the calculation. We are adding the 
coronal approximation to treat such lines. 

There are many molecular lines in the observed spectrum. We are adding 
molecular lines to our line list, 

Eor many lines lilce Fe II we use an approximate scattering source function in 
the chromosphere, but emission cores still appear in the calculation, contrary to 
observation. Wo are developing a more accurate procedure for dealing with these 
lines. 
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